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Abstract 

This  paper  analyzes  and  compares  the  performance  of  two  prominent  explicit  guidance  algorithms:  the  iterative  guidance  mode  and 
the  powered  explicit  guidance.  We  performed  a  series  of  numerical  simulations  of  a  space  launch  vehicle  model  for  both  nominal  and 
off-nominal  conditions.  One  of  our  findings  is  that  if  we  take  into  account  the  originally  ignored  higher-order  terms  from  the  guidance 
parameters  of  the  iterative  guidance  mode  for  a  long-range  flight,  the  guidance  performance  can  be  enhanced  to  a  level  comparable  to 
that  of  the  powered  explicit  guidance.  These  higher-order  terms  can  be  included  by  employing  an  iterative  predictor-corrector  method 
like  the  powered  explicit  guidance.  Also  we  proposed  a  remedy  of  preventing  relatively  earlier  divergence  of  the  guidance  commands  with 
the  predictor-corrector  iteration  than  that  of  the  linear  differential  corrector  approach  by  making  a  better  initial  guess. 
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Keywords:  Explicit  guidance  law;  Space  launch  vehicle;  Thrust  acceleration  integrals;  Iterative  methods;  Guidance  performance  analysis 


1.  Introduction 

Most  of  the  modern  multi-staged  space  launch  vehicles 
adopt  closed-loop  guidance  laws  for  the  exoatmospheric 
trajectory,  although  open-loop  ascent  guidance  is  effective 
enough  for  the  first  stage  maneuver  against  the  atmo¬ 
spheric  loads.  One  of  the  main  virtues  of  a  closed-loop 
guidance  program  is  that  the  steering  vectors  and  engine 
cutoff  times  can  be  determined  to  achieve  the  target  with 
minimum  fuel  consumption.  For  this  purpose,  various 
explicit  guidance  schemes  based  on  linear  tangent  laws 
have  been  developed  since  the  space  era  of  1960s  and  their 
inherent  mission  flexibility  has  allowed  them  to  be  applied 
even  in  today’s  space  launch  vehicles  (Shrivastava  et  al., 
1986).  To  handle  more  complicated  mission  requirements 
for  future  space  transportation,  guidance  algorithms  with 
higher  degree  of  onboard  adaptability  have  been  being 
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developed  (Hardtla  et  al.,  1987;  Jessick  and  Knobbs, 
1992),  including  those  based  on  onboard  trajectory  optimi¬ 
zation  with  advances  in  computer  processors  (Skalecki  and 
Martin,  1993;  Jutty  et  al.,  2000).  Other  research  efforts,  also 
principally  indebted  to  much  improved  onboard  computer 
capabilities,  have  been  devoted  to  eliminate  simplifying 
assumptions  previously  required  to  derive  the  closed-form 
solutions  of  guidance  problems  (Delporte  and  Sauvient, 
1992). 

It  could  be  said  that  the  iterative  guidance  mode  (IGM) 
developed  for  guidance  of  the  Saturn  launch  vehicles 
(Chandler  and  Smith,  1967)  and  the  powered  explicit 
guidance  (PEG)  algorithm  developed  during  the  Space 
Shuttle  Program  with  various  missions  (McHenry  et  al., 
1979;  Schleich,  1982)  are  among  the  most  distinguished 
and  flight-proven  guidance  algorithms  ever  devised. 
Although  the  performance  of  these  algorithms  has  been 
demonstrated  by  hundreds  of  successful  flights,  it  is  rather 
hard  to  find  open  literatures  analyzing  the  effects  of  the 
approximations  used  in  their  derivation  or  comparing  the 
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performance  of  these  two  algorithms.  Song  et  al.  (2011) 
proposed  a  method  to  reduce  the  prediction  errors  of 
IGM  by  taking  into  account  the  effect  of  approximate 
formula  of  the  angle-to-go  prediction  to  orbit  injection 
accuracy.  A  higher-order  form  of  PEG  was  proposed  by 
Sinha  and  Shrivastava  (1990)  to  consider  the  high  maneu¬ 
vers  of  Indian  launch  vehicles,  and  they  showed  that  its 
performance  was  superior  to  that  of  the  original  form. 
Vittal  and  Bhat  (1991)  and  Delporte  and  Sauvient  (1992) 
formulated  the  thrust  acceleration  integrals  with  constant 
thrust  analytically  for  the  vector  form  of  the  linear 
tangent  law,  which  had  been  originally  approximated  by 
Taylor  series  values  (McHenry  et  al.,  1979;  Sinha  and 
Shrivastava,  1990). 

This  study  analyzes  and  compares  the  performances  of 
IGM  and  PEG  while  seeking  methods  to  improve  the 
performance  of  these  guidance  schemes  in  their  original 
forms.  We  investigated  the  effect  of  higher-order  nonlinear 
terms  previously  ignored  in  IGM  to  derive  the  approximate 
linear  equations  of  guidance  parameters  by  employing 
similar  approaches  by  Sinha  and  Shrivastava  (1990)  with 
PEG  algorithm.  That  is,  higher-order  terms  are  shown  to 
be  included  in  almost  the  same  way  with  PEG  algorithm 
if  an  iterative  predictor-corrector  method  is  adopted. 
There  are  some  references  commenting  that  PEG  can  be 
considered  as  a  kind  of  a  vector  form  of  IGM  and  the 
performance  of  these  two  algorithms  is  nearly  identical 
for  vacuum  flight  phase  (Jaggers,  1977;  Hanson  et  al., 
1995).  However,  simulation  results  show  that  the  perfor¬ 
mance  of  IGM  can  be  more  deteriorated  than  that  of 
PEG,  especially  in  terms  of  slower  convergence  in  the 
time-to-go  prediction.  It  is  shown  that  the  performance 
of  IGM  can  be  improved  by  introducing  previously 
ignored  higher-order  terms,  which  again  makes  the  linear 
angle  command  assumption  of  the  algorithm  more  valid. 

The  guidance  parameters  of  PEG  are  computed  using 
the  iterative  predictor-corrector  approach  starting  with 
an  initial  guess  at  every  major  guidance  computation  cycle. 
By  comparing  with  the  linear  differential  correction  method 
using  first  derivatives  to  all  guidance  parameters  (Delporte 
and  Sauvient,  1992),  it  is  shown  that  the  initial  guess  of 
the  predictor-corrector  approach  is  more  erroneous  as 
the  time-to-go  approaches  0,  which  may  lead  to  more  iter¬ 
ations.  It  is  well  known  that,  for  iterative  algorithms,  the 
first  guess  close  to  the  answer  significantly  improves  the 
convergence  properties  (Luenberger  and  Ye,  2008). 
The  same  applies  to  the  initial  guess  of  the  velocity-to-go, 
the  running  variable  of  the  iteration  process  in  PEG.  In  this 
study,  we  showed  that  the  convergence  of  the  iteration 
process  in  PEG  can  be  enhanced  by  improving  the  initial 
guess  of  the  velocity-to-go. 

Although  our  approaches  here  may  not  claim  signifi¬ 
cant  performance  improvement  on  these  well-developed 
schemes  in  their  original  form,  we  proved  that  some  defi¬ 
ciencies  can  be  overcome  without  much  more  complexities 
such  as  onboard  optimization  methods  which  are  currently 
one  of  the  major  research  topics  in  the  space  launch  vehicle 


guidance  area.  Except  for  this  introduction,  this  paper  is 
organized  into  four  sections.  The  formulae  of  higher-order 
form  of  IGM  are  elaborated  first  and  then  described  is  a 
comparison  between  the  predictor-corrector  approach  of 
PEG  and  the  linear  differential  corrector  method.  Follow¬ 
ing  the  simulation  results  for  the  proposed  approaches, 
the  conclusions  of  this  work  are  summarized. 

2.  Explicit  guidance  law  with  higher-order  nonlinear  terms 

Most  of  the  explicit  guidance  laws  have  a  functional 
form  based  on  an  optimal  trajectory  minimizing  the  usage 
of  propellants.  Derivation  of  IGM  is  also  based  on  the  fact 
that  the  optimal  thrust  direction  for  the  flat-earth  and 
uniform  gravity  model  is  approximately  given  by  a  linear 
angle  form  (Chandler  and  Smith,  1967).  The  pitch  and 
yaw  attitude  commands  are  assumed  to  be  linear  relative 
to  the  guidance  reference  system  G  as  follows,  with  G  being 
defined  in  the  target  orbit  plane  whose  origin  is  located  at 


the  center  of  the  Earth. 

9g  =  9Gv  — 

(oGp  — 

(1) 

-<5- 

II 

(<A Gr  -  bat ) 

(2) 

where  current  time  is  set  to  0,  and  t  denotes  the  time  vary¬ 
ing  from  0  to  time-to-go,  tgo,  the  remaining  time  before 
achieving  the  target  injection  point.  The  x-axis  of  the  guid¬ 
ance  coordinate  system  G  goes  toward  the  predicted  orbit 
injection  point  from  the  origin,  the  y-axis  hes  in  the  normal 
direction  of  the  target  orbit  plane,  and  the  --axis  completes 
the  right-handed  coordinate  system.  Among  the  guidance 
parameters,  dGv,0Gp ,  0G,  ^Gv,^G  ,  and  i//G  are  related  to  the 
commanded  attitudes  (thrust  direction)  and  tgo  determines 
the  time  of  the  engine  shutdown.  Finally  range-angle-to-go, 
tgo,  an  angle  between  current  position  and  predicted  target 
injection  point,  should  be  computed  to  define  the  G  frame. 
These  parameters  are  calculated  to  satisfy  the  position  and 
velocity  constraints  of  a  target  orbit  by  analytic  integration 
of  the  following  point-mass  equations  of  motion  applicable 
to  a  typical  space  launch  vehicle  in  a  vacuum. 


i  sin  0G  cos  i j/G  +  g G 

(3) 

i  sin  i j/G  +  gGy 

(4) 

!  COS  0G  COS  1 l/G  +  gf 

(5) 

where  the  first  terms  represent  the  acceleration  by  thrust 
and  the  second  terms  by  gravity,  which  are  integrated 
independently. 

The  derivation  procedure  of  the  guidance  parameters 
for  IGM  in  its  original  (Chandler  and  Smith,  1967)  and  a 
modified  form  to  consider  previously  ignored  higher-order 
terms,  is  as  follows.  First,  dGv  and  i j/Gv  are  derived  to  satisfy 
only  the  terminal  velocity  constraints  under  the  assumption 
of  invariant  thrust  direction.  Integration  of  Eqs.  (3)— (5) 
during  the  remaining  burn  time  leads  to 
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xf  =  xG  +  sin  0Gv  cos  t i//Gp  J  adt  +  vGrav* 

(6) 

yG  =  yG  +  sin  i l/Gv  J  adt  +  vGrav^ 

(7) 

zG  =  zG  +  cos  0G„  cos  \jjGr  j  adt  +  vGrar  . 

(8) 

Then  the  guidance  parameters  are  obtained  as  follows. 


9Gv  =  tan  1  (  -r-^  ,  i j/G  =  tan  1  -  = 

\Jav?  +  A  vf 


(9) 


where 

Av°  =  (10) 

Here  tgo  is  derived  from  the  following  relationship. 

(11) 


Then  the  remaining  attitude  guidance  parameters 
9G ,  0G,  1 1>G  and  il/G  are  derived  to  satisfy  the  terminal  posi¬ 
tion  constraints  in  the  directions  of  x  and  y  axes.  To  derive 
analytic  expressions  for  the  thrust  integrals,  ( 9Gp  —  0Gt )  and 
(•Ag  -  <AgO  are  assumed  to  be  small  enough  so  that  the  sine 
and  cosine  functions  of  t  can  be  substituted  by  series 
approximation  as  follows. 


sin(0qp  —  SGt)  ( 0Gp  —  0Gt )  - 


;p-eGt)3 


(12) 


sin  (<Agp  -  •Aci)  ^  G"  „  ^  +  ■ 

/  .  \  (V Gp  —  &Gt)  (*Agp  —  ^Gt) 

cos  (VGp  -  j,Gt)  « I-  V  2,  y  V  4,  y 


(13) 

(14) 


(15) 


The  original  IGM  algorithm  considers  only  the  first 
underlined  terms  in  Eqs.  ( 1 2)— ( 1 5)  and  such  ignorance  of 
higher-order  terms  is  valid  only  when  time-to-go  is  suffi¬ 
ciently  small.  From  these  disregarded  terms,  we  can  infer 
that,  for  long  range  cases,  the  thrust  integrals  may  not  be 
accurate  and  even  deteriorate  the  guidance  performance. 
According  to  the  simplification  adopted  in  the  original 
IGM,  the  x-axis  component  of  thrust  acceleration  in 
Eq.  (3)  can  be  obtained  as  follows. 

sin  0G  =  sin  ( 'oGv  —  (oGp  —  0g/)) 

«  sin  dGv  -  (oGp  -  9Gtj  cos  6Gv  (16) 


cosine  =  cos  (VGc  -  (t//Gp  -  <Ag*)) 

«  COS  i l/Gv  +  (Vgp  -  sin  i//Gv  (17) 

sin 0G cos  \j/G  =  sm9Gv  cosi/tGt  -  (oGp  -  9Gt^j  cosdGv  cost \j/Gc 
+  (j/Gp  -  iAc?)  sin0Go  sin  i /rGs 
-  (9Gp-  ('A Gp  -  'Ac?)  COS  0G„  sin  I l/Gv  (18) 

It  can  be  noticed  that  the  original  IGM  does  not  con¬ 
sider  all  the  second-order  terms  in  an  exact  manner.  From 
Eqs.  ( 1 6)— ( 18)  it  may  appear  that  all  terms  less  than  or 
equal  to  the  second-order  terms  of  t  are  considered,  how¬ 
ever,  if  the  terms  with  double  underlines  in  Eqs.  (13)  and 
(15)  are  not  neglected  from  the  first  step  of  the  algorithm 
derivation,  then  as  in  Eq.  ( 19),  another  second-order  term 
should  show  which  may  have  equal  importance  to  other 
second-order  terms  in  Eq.  (18).  Likewise  the  third-  and 
more  higher-order  terms  also  need  to  be  included  as  shown 
in  Eq.  (19)  for  accuracy  enhancement. 

sin0G  cost/iG  =  sin  Do,  coSI/rG,  -  (Bo,  -  0ot )  cos0G,  cosi/«G.  +  (j,Gp  -  <m)  sin0G,  sin ^ 
-  (0O,  -  6oi)  (> K,  ~  Ac')  cos0G.  sini/rG- 
-\{6o,-0at)  sm0G„cosi/iG,  “  ^  (Ao,  “  Ag<)  sin0G,cosi/ro< 

+~  (eGp  -  0Gr)  2  (i/r Gr  -  j/at) 2  sin0G,  cosi •  •  •  (19) 

Using  Eq.  (19),  the  x-axis  velocity  component  is 
obtained  by  integration  of  Eq.  (3)  during  tgo  and  additional 
integration  gives  the  x-axis  position  component  as  follows. 

ief + j  a  sin  0G  cos  icdt+v^ 


=if +Isin0G,  cosi/<Gi  -  (z0G,  -70e)  cos 0O,  cosi^o,  +  (LfCt  -jjie)  sin0G,  sin^G> 


-Be,  (LiGp  -Ji^g)  cos  0g, 

sin^a+0G(^Gf-^G)c< 

>s0G,sini/ia 

-i  (ifgg,  -  UB0,  Be  +H82c 

)  sin0G,cosi/ift  -\[lVgp sini 

9Gb  cosi !/Gb 

+^sin0G,  cosi/rG,{£0^,iAG, 

-2/(02G,^AG  +  0Gp0G^; 

(20) 

(=xG0+x%0  +  rfa^De, 

jos  \jiGdsdt  +  p^rav 

=x^+xGtgo+SsmBG,  cos 

-  (seGr  -  QBg)  cos0g,  cosi/iG, 

+(^<i-Assr)«fa 

Be,<m*o, 

-Be,  (S<PC,  ~  fiAs)  cos0G,  sin +0G  (&,Gp  -P^G)  cos 0C 

-lisfto-ZQBcA+PPa)* 

mBa,cos*0p-\(stfh,-2Qtl 

3>g+pAo)s inflao 

081 l>G. 

+isin0G,  cosi Ag„  {Stfe,  l/fG,  -i 

tQ^K^+B^e^e,)- 

+pI«, 

(21) 

where  L,  J,  H,  S,  Q,  and  P  are  defined  as  follows. 

L  =  f‘os°  adt  S  =  J‘os°  J‘o  adsdt 

J  =  f‘os°  atdt  Q  =  j‘os°  j‘o  asdsdt  (22) 

H  =  J‘os°  at2dt  P  =  J‘gs°  f‘o  as2dsdt 

For  analytic  integration  of  Eq.  (22)  over  the  burning 
time,  constant  thrust  or  constant  acceleration  models  are 
usually  employed.  In  a  similar  manner,  for  the  y-axis  veloc¬ 
ity  and  position  components, 


Please  cite  this  article  in  press  as:  Song,  E.-J.,  et  al.  A  comparison  of  iterative  explicit  guidance  algorithms  for  space  launch  vehicles.  J.  Adv.  Space 
Res.  (2014),  http://dx.doi.Org/10.1016/j.asr.2014.09.025 


ARTICLE  IN  PRESS 


4  E.-J.  Song  et  al.  I  Advances  in  Space  Research  xxx  (2014)  xxx—xxx 


yGf=yG0+  / 

Jo 

a  sin  i jjGdt  +  vgrar> 

=  yGa  +L  sin 

'Ac,,  -  (Lll'Gp  ~  J^(j) 

COS  1 j/Gv 

-jK 

-2JiAg,<K?  +  #iAg)s 

iin  'Ac,  •  • 

•  +  Vtavy  (23) 

yGf=yG0+y%° 

,  +  J  J  a  sin  \j/Gdsdt  +  pgrai 

=  yG+y%o 

+  S  sin  i j/Gti  -  (S\l/Gp 

-  e<Ac) 

COS  lAct 

-2Q'I'Gp'I'g+P'I'1J)  i 

iiniA  c,  •' 

■+P%aVy  (24) 

Many  real  mission  designers  have  taken  the  great  possi¬ 
ble  care  to  avoid  large  yaw  maneuvers  to  change  the  orien¬ 
tation  of  the  orbital  plane  since  these  kinds  of  maneuvers 
are  not  economical  in  terms  of  propellant  consumption. 
In  that  case  the  higher-order  terms  of  (i AG  -  \j/Gt)  in  the 
yaw  motion  can  be  neglected.  Then  the  four  equations, 
Eqs.  (20),  (21),  (23)  and  (24),  are  used  to  compute  the  four 
unknowns  dGp,  0G,  t j/G  ,  and  i//G,  and  the  velocity  equations 
of  Eqs.  (20)  and  (23)  come  to  the  same  equations  as  Eqs. 
(6)  and  (7),  respectively,  from  the  following  equations. 

—  (l8gp  —J0g\  cos 0g,  cosi AG>  +  (liI>gp  —J^g)  sin0G,  siniAG> 

-eGp  (l*Gp  -aa  g)  COS0G,  simAG-  +ea  (aag, -a  h)  cos  ea,  sin*G, 

~  (l^,  -  Z I0Gp  i,G+Hi?G)  sin0G,  cosi/rG,  -A  (l^  -  ZJ^Gp i/rG  sin0G,  cosi AG, 

-4sinflG,cosiAGi{i^I,iAG,-2J^piAG,iAG+0G,0Gi^,)-'  =  O  (25) 

-(l'I'gp-j'I'g)cosiI/Gv- -  (liI/2Gp  -2J\I/Gpil/G+ sin i/f Gv  ■  ■  ■  =  0 

(26) 

For  the  velocity  in  z-axis  direction, 


=zf  +Zcos0G,  cosi Ag,  +  (l6Gp  -J8a)  sin0G„  cosi/rG,  +  (ii/«Gj,  -  AAG)  cos0G,  siniAG- 

+  0Gp  (Li//Gp  — A//g)  sin  0Gl  sin  ij/fg  —  0G  ( J<l/Gp  —  Hi) /G)  sin  0Gl  sin  i//Gt 

-  i  (Z0^  -  ZJeGp  0G + Hi 0* )  cos  0g.  cos  I AG„ 

-  i  (il/rG,  -  ZJ>l>Gp  iAg  +  ffi>G)  cos  0C,  COS  I Ac, 

+  icosflG, COS lAc, {-£00,^0, -Zj(92Gpil/Gpil/a+9Gp()o<l'lr>)  ■  ■  +  (27) 

For  the  original  IGM,  approximate  linear  equations  are 
derived  for  the  four  parameters  9Gp,6G,\l/G,  and  t j/G  by 
ignoring  the  higher  order  terms.  In  that  case,  Eq.  (27)  also 
becomes  the  same  as  Eq.  (8).  The  z-axis  position  equation 
can  be  derived  in  a  similar  manner,  where  the  constraint 
value  is  0  from  the  definition  of  the  G  frame.  The  guidance 
parameter,  cg0  can  be  computed  from  the  z-axis  position 
equation,  however,  Chandler  and  Smith  (1967)  introduced 
a  mission-dependent  constant  Kg  to  calculate  £,go  approxi¬ 
mately  to  avoid  solving  the  z-axis  position  equation,  as 
follows. 


Zgo  =  l(vtgo  +  S  +  Kgt2go)  (28) 

where  Kg  is  calculated  before  each  flight  through  the 
computer  simulations  of  the  flight  trajectories. 

The  guidance  algorithm  is  a  procedure  to  solve  the 
equations  for  the  guidance  parameters.  The  unknown 
parameters  are  0G:: ,  9Gp ,  9G.  ^tGr ,  ijjG  ,  \j/G-  tg0,  and  tgo  (8 
unknowns),  and  they  must  satisfy  the  target  position  and 
velocity  constraints  (6  equations)  and  two  additional  veloc¬ 
ity  constraints  (2  equations)  to  simplify  the  procedure  of 
finding  solutions  by  formulating  the  approximate  linear 
equations  (Chandler  and  Smith,  1967).  In  that  case,  with 
given  tgo,  the  iteration  is  performed  only  for  tgo  and  Qgo. 
Then  the  guidance  equations  consisting  of  only  the 
first-order  and  some  selected  second-order  terms  of  the 
guidance  parameters  related  to  maneuvers  can  be  solved 
analytically. 

On  the  other  hand,  in  order  to  consider  the  higher-order 
nonlinear  terms,  the  iteration  must  include  the  maneuver- 
related  guidance  parameters,  and  thus  is  accompanied  by 
a  little  more  computational  complexity.  To  solve  nonlinear 
equations  with  higher-order  terms,  either  the  predictor- 
corrector  approach  or  the  linear  differential  method  can 
be  employed.  Between  these  two  of  the  most  prevailing 
solvers,  we  chose  the  former  approach  here  as  in  the  origi¬ 
nal  PEG.  A  predictor-corrector  algorithm  consists  of  two 
separate  steps  to  find  a  solution:  the  prediction  step  to 
calculate  an  approximation  of  the  desired  quantity,  and 
then  the  correction  step  to  refine  the  initial  approximation 
using  a  different  method.  Likewise,  in  our  approach,  the 
initial  guidance  parameters  are  computed  by  solving  the 
approximate  linear  equations  as  in  the  case  of  the  original 
IGM  and  then  the  higher-order  terms  are  computed  using 
the  approximate  solution.  In  the  next  iteration  step,  the 
higher-order  terms  ignored  are  added  to  the  linear 
equations  as  constant  bias  terms.  In  this  way,  the  solutions 
can  be  improved  as  the  iteration  proceeds.  A  major 
advantage  of  the  original  IGM  is  that  it  does  not  require 
massive  iterations,  at  the  same  time  it  causes  limited  upper 
bounds  for  the  applicable  region  (Karacsony  and  Cole, 
1970).  By  including  the  higher-order  terms  as  proposed  in 
this  work,  the  applicability  of  IGM  can  be  easily  extended 
and  the  corresponding  computational  cost  is  expected  to  be 
similar  to  that  of  PEG. 


3.  Explicit  guidance  law  using  predictor-corrector  algorithm 

The  guidance  parameters  of  PEG  are  slightly  different 
from  those  of  IGM,  since  the  thrust  direction  is  given  by 
a  simplified  vector  form  of  the  linear  tangent  steering  law 
as  shown  in  Eq.  (29),  which  is  derived  from  an  ascent  tra¬ 
jectory  solution  with  minimum  fuel  consumption 
(McHenry  et  al.,  1979). 
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UF 


X„  +  A(t  —  tx) 


~  AV{1  —  2  'X2^  —  ft)2}  +  X(t  —  tx ) 


where 


iX)i  =  l)  i  i=o 


Convergence  can  be  determined  based  on  the  updated 
time-to-go  not  changing  significantly  from  that  of  the 
previous  iteration  step.  . 

The  foregoing  seven  guidance  parameters  Av,l,  and  tgo 

(29)  of  PEG  can  also  be  computed  using  the  linear  differential 
correction  method  as  proposed  by  Delporte  and  Sauvient 
(1992).  The  nonlinear  relationship  between  the  guidance 
parameters  and  the  desired  cutoff  constraints  can  be  repre- 

(30)  sented  as  follows. 


As  briefly  mentioned  in  the  previous  section,  the  predic¬ 
tor-corrector  computational  sequence  is  adopted  to  solve 
the  guidance  parameters  satisfying  the  cutoff  constraints 
in  PEG  algorithm.  The  seven  unknown  parameters  AV,A, 
and  tgo  are  computed  in  iterative  fashion  using  v go  as  a 
running  variable.  The  velocity-to-go  of  each  major  loop 
is  calculated  to  be  diminished  by  the  velocity  change  At; 
from  the  vehicle  velocity  of  the  previous  iteration  cycle. 

vgo(n)  =  vgo(n  -  1)  -  AtJ  (31) 

With  given  vgo,  the  prediction  process  starts  with  evalu¬ 
ation  of  tgo  and  the  thrust  integrals  defined  in  Eq.  (22). 
Then  the  first  and  second  integrations  of  the  point-mass 
equations  of  motion  using  the  thrust  direction  form  in 
Eq.  (29)  give  the  following  linear  relationships  among  the 
guidance  parameters  and  cutoff  conditions. 


(32) 

j_rgo-sl 
( Q-Stx ) 

(33) 

j 

h~L 

(34) 

where 

rgo  =  rd-r-  vtgo  -  rgrap  +  rbias 

(35) 

The  bias  term,  rbias  =  rgo  —  fthmst  ats  \  A2  (P  —  2 Qtx  +  Sify 
Av,  is  computed  using  the  guidance  parameters  obtained 
from  the  previous  iteration  step.  In  this  way,  the  second- 
and  higher-order  nonlinear  terms  can  be  included  as  a  bias 
term  and  considered  in  the  prediction  of  the  cutoff  velocity 
and  position. 


Vp  =  V  +  Vthrust  +  Vgrav 

(36) 

rp  =  r  +  vtgo  +  f thrust  +  rgmv 

(37) 

where 

vthrus,^Av{L~A\H-Jtx)  J 

(38) 

thrust  «  -  \'A2{P  -  2  Qtx  +  Stl)  | 

+  (Q  -  Stx)A. 

(39) 

Next  in  the  correction  process 

,  the  desired 

cutoff 

condition  and  the  velocity-to-go  ; 

are  computed 

using 

the  predicted  states  and  the  target  orbit  information. 


Y=f(U)  (40) 


where 


1 

0 


The  former  five  constraints  from  Eq.  (41)  represent  the 
desired  position  and  velocity,  which  are  usually  given  by 
target  orbital  parameters,  and  the  latter  two  constraints 
are  derived  from  Eq.  (30).  The  guidance  parameters  U 
are  updated  until  the  targeting  error  AT  becomes  close  to  0. 


ir=m&u 


(42) 

(43) 


The  difference  between  the  predicted  end  points  and  the 
targets,  AT,  can  be  computed  by  analytical  integration  of 
the  flight  trajectory  in  the  same  way  as  for  PEG  starting 
with  an  initial  guess  of  U  given  from  the  previous  guidance 
cycle.  The  sensitivity  matrix,  can  also  be  evaluated 
analytically  from  the  nonlinear  relationship  between  U 
and  T.  Detailed  description  of  the  derivatives  of  the 
constraints  with  respect  to  the  guidance  parameters  are 
presented  in  Delporte  and  Sauvient  (1992). 

Song  et  al.  (2013)  compared  the  predictor-corrector 
approach  of  PEG  with  the  linear  differential  corrector 
method  and  showed  that  the  guidance  parameters  obtained 
from  both  methods  became  very  similar  if  the  thrust 
and  gravity  acceleration  integrals  needed  in  trajectory 
integrations  were  computed  in  the  same  way.  However 
their  iterative  algorithms  described  above  imply  that, 
regardless  of  the  similarity  between  the  converged  solu¬ 
tions,  intermediate  solutions  from  the  iteration  processes 
can  be  slightly  different.  In  PEG,  the  decrement  of  the  inde¬ 
pendent  variable  vgo  from  previous  guidance  cycle  is  not 
performed  in  an  exact  manner,  that  is,  only  sensed  change 
in  velocity  is  considered  with  the  start  of  current  iteration 
cycle  as  in  Eq.  (31).  Moreover,  inaccuracy  of  the  initial 
guess  can  induce  more  iteration  and  even  early  divergence 
of  the  guidance  commands,  in  other  words,  guidance  insta¬ 
bility.  Most  of  the  explicit  guidance  schemes  for  launch 
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vehicles  including  PEG  have  a  general  stability  problem, 
especially  when  the  vehicle  is  close  to  the  injection  point. 
For  PEG,  the  denominators  of  Eqs.  (32)  and  (33)  approach 
0  as  tgo  reduces  to  0,  which  explains  the  intrinsic  instability 
of  the  scheme.  Bittner  (1976)  and  Vittal  and  Bhat  (1993) 
have  tried  to  solve  the  problem  by  formulating  a  guidance 
law  using  a  quadratic  injection  criterion  rather  than  mak¬ 
ing  injection  errors  to  be  exactly  zero.  These  previously 
suggested  approaches  have  focused  on  eliminating  undesir¬ 
able  forms  of  the  guidance  parameters  equations  which 
make  the  system  unstable;  we  treated  this  problem  in 
another  way  by  taking  into  account  the  effect  of  the  itera¬ 
tive  algorithm  for  PEG  to  guidance  stability.  Although  the 
complete  avoidance  of  the  instability  is  not  feasible  with 
the  identical  equations  for  the  guidance  parameters,  it 
can  be  delayed  by  selecting  initial  guess  more  elaborately. 
Luenberger  and  Ye  (2008)  have  described  that  the  solutions 
of  many  iterative  algorithms  for  nonlinear  problems  have 
sensitivity  on  the  starting  points,  and  we  noticed  that  this 
kind  of  effect  may  become  noticeable  as  the  time-to-go 
approaches  0.  If  changes  in  the  gravity  integrals  are 
additionally  considered,  Eq.  (31)  turns  into  Eq.  (44),  and 
this  slight  modification  can  correct  the  inaccuracy  caused 
by  the  predictor-corrector  approach. 

vgo  (n)  =  vgo  (n  —  1 )  —  Av  —  Avgrav  (44) 

where 

Atgrav  =  Vgrav  («)  -  Vgrav  («  -  1 )  (45) 

The  gravity  integral  terms  can  be  easily  calculated  using 
navigated  states  information  in  the  current  step  and  the 
predicted  final  states  obtained  from  the  previous  cycle. 


For  example,  when  using  the  method  proposed  by 
Delporte  and  Sauvient  (1992)  and  Jaggers  (1977),  the  grav¬ 
ity  integral  is  given  by 

Vgrav  =  y  (1/  +  g0)  +  ^  (g„  ~  g/j  (46) 

which  is  derived  from  the  gravity  field  approximated  by  a 
third-degree  polynomial  function  such  as 

g  =  E<  (47) 

The  coefficients  a0,  a\,  a2,  and  a3  are  evaluated  to  satisfy 
the  gravity  and  its  first  derivative  at  current  position  and 
the  predicted  final  position,  g0,g/,g0,  and  gf.  More  details 
about  this  approach  are  provided  in  Delporte  and  Sauvient 
(1992)  and  Jaggers  (1977).  Now  vgmv(n)  in  Eq.  (45)  to 
obtain  Avgrm  is  evaluated  as  follows. 

+  (aw - !/(» - 1))  (48) 

where  AT  denotes  the  elapsed  time  from  the  previous  cycle. 
Then  the  convergence  characteristics  can  be  improved  to  a 
similar  level  with  the  linear  differential  corrector,  where  AT 
can  be  considered  exactly  in  the  variable  tgo  as  shown  in 
Eq.  (49). 

tgo{n)  =  tgo(n  —  1)  —  AT  (49) 


Flight  Time  (sec) 


Fig.  1.  Pitch  attitude  command  during  ascent  flight  generated  by  different  guidance  schemes. 
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Fig.  2.  Pitch  attitude  command  close  to  the  orbit  injection  point  generated  by  different  guidance  schemes. 


Other  guidance  parameters  lv  and  1  in  Eq.  (41)  change 
slowly  during  AT,  which  means  their  initial  guesses  using 
the  results  from  the  previous  iteration  cycle  are  so  close 
to  the  solutions  from  the  current  cycle  and  they  do  not  trig¬ 
ger  any  convergence  problems.  Comparison  between  the 
original  approach  and  proposed  modification  is  presented 
in  the  next  section. 


Table  1 

Nominal  trajectory  data  for  the  3-staged  launch  vehicle  model. 


First  stage 

Second  stage 

Third  stage 

Ignition 

time  (sec) 

0 

12S 

276 

Burn  tir 

ne  (sec) 

126 

144 

502 

Flight  Time  (sec) 


Fig.  3.  Yaw  attitude  command  in  stage  3  generated  by  different  guidance  schemes. 
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Fig.  4.  Prediction  error  of  tgo  with  the  nominal  and  off-nominal  conditions.  ‘IGM’  shows  noticeably  larger  errors  than  those  of  other  schemes. 
‘IGM  +  H.O.T.’  fixes  the  problem  and  shows  similar  performance  with  ‘PEG’. 


4.  Guidance  performance  evaluation 

We  performed  a  series  of  3-DOF  (degree  of  freedom) 
computer  simulations  here  for  an  ascent  trajectory  of  a  3- 
staged  launch  vehicle  model  with  following  five  closed-loop 
guidance  algorithms  applied  during  the  exoatmospheric 
phase  of  the  second  and  third  stage  flights. 

(1)  IGM  (original) 

(2)  IGM  +  H.O.T.  (higher-order  terms) 

(3)  PEG  (original) 


(4)  PEG  +  M.I.G.  (modified  initial  guess) 

(5)  PEG  +  L.D.C.  (linear  differential  corrector) 

The  closed-loop  guidance  operation  is  initiated  after 
about  10  (sec)  from  the  ignition  of  each  upper  stage.  Before 
then  open-loop  commands  are  employed  just  like  for  the 
stage  1  flight,  where  the  vehicle  ascends  against  the  dense 
atmosphere.  It  is  usual  to  adopt  an  open-loop  logic  for  this 
phase  for  safety  and  reliability  reasons.  Target  orbits  of  the 
stage  2  and  stage  3  were  set  to  be  their  nominal  final  con¬ 
ditions  obtained  by  trajectory  optimization  (‘Optimal’  in 


Table  2 

Stage  3  shutdown  sensitivities  to  the  nominal  condition  with  different  guidance  algorithms. 


Guidance  algorithm 

R(Perigee)  (m) 

R( Apogee)  (m) 

Inclination  (deg) 

Longitude  of  ascending  node  (deg) 

vgo  <  150  (m/s) 

IGM 

143.08 

-306.61 

-0.00025 

0.00005 

IGM  +  H.O.T. 

99.38 

-297.31 

0.00072 

-0.00012 

PEG 

182.28 

-287.85 

-0.00010 

0.00002 

PEG  +  M.I.G. 

182.25 

-288.74 

-0.00010 

0.00002 

PEG  +  L.D.C. 

182.27 

-287.87 

-0.00010 

0.00002 

Vgo  <  500  (m/s) 

IGM 

-1488.84 

1347.83 

0.00259 

-0.00047 

IGM  +  H.O.T. 

-100.32 

-485.51 

0.00162 

-0.00028 

PEG 

-238.50 

-392.96 

0.00118 

-0.00021 

PEG  +  M.I.G. 

-238.45 

-393.00 

0.00118 

-0.00021 

PEG  +  L.D.C. 

-238.73 

-393.02 

0.00118 

-0.00021 

vgo  <  750  (m/s) 

IGM 

-4201.31 

3403.95 

0.00865 

-0.00161 

IGM  +  H.O.T. 

-629.79 

-347.72 

0.00387 

-0.00070 

PEG 

-1228.23 

-283.23 

0.00580 

-0.00106 

PEG  +  M.I.G. 

-1228.13 

-283.31 

0.00119 

-0.00021 

PEG  +  L.D.C. 

-1228.06 

-283.18 

0.00580 

-0.00106 
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Fig.  5.  Scatter  plot  of  insertion  perigee  error  versus  apogee  error  depending  or 
Earlier  guidance  stop  with  ‘IGM’  results  in  larger  insertion  errors  than  those  of  < 
similar  results. 


i)  vgo  s/  150(m/s)  and  (b)  vgo  sc  750(r 
I.I.G.’,  and  ‘PEG  +  L.D.C.’  provide 


Figs.  1-3).  After  the  burnout  of  the  third  stage  engine,  the 
payload  is  assumed  to  be  injected  into  the  target  orbit  spec¬ 
ified  as  a  700  km-high  near-circular  orbit. 

To  compare  above  guidance  algorithms,  here  we  adopt 
the  same  gravity  integrals  prediction  method  that 
Delporte  and  Sauvient  (1992)  and  Jaggers  (1977)  have 
employed  in  their  work.  As  briefly  described  in  Section  3, 
the  gravity  field  is  modeled  with  a  third-degree  polynomial 
function  of  time.  To  ascertain  the  effect  of  the  higher-order 
terms  described  in  Section  2,  the  flight  duration  of  the  third 


stage  is  set  to  be  much  longer  than  that  of  the  second  stage 
for  the  nominal  trajectory,  as  shown  in  Table  1.  We  evalu¬ 
ated  the  guidance  performance  for  off-nominal  conditions 
as  well  as  for  the  nominal  condition  to  validate  the  sug¬ 
gested  methods  by  conducting  the  Monte  Carlo  analysis. 
The  Monte  Carlo  analysis  is  a  practically  useful  tool  to 
evaluate  the  effects  of  complicated  combinations  of  multi¬ 
ple  variables  and  broadly  employed  for  launch  vehicle 
design  and  performance  analysis  (Bell  et  al.,  1993; 
Hanson  and  Hall,  2008).  For  each  guidance  algorithm, 
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Fig.  6.  Yaw  attitude  command  in  the  guidance  frame  of  IGM  with  the  nominal  and  off-nominal  conditions.  The  linearity  assumption  is  more  valid  with 
‘IGM  +  H.O.T’. 


the  Monte  Carlo  analysis  consisting  of  500  individual  sim¬ 
ulations  was  performed  with  simultaneously  varying  error 
sources  modeling  the  main  performance  uncertainties  in 
thrust,  specific  impulse,  dry  and  propellant  weights,  and 
so  on  for  each  stage.  Here  each  error  source  was  modeled 
as  a  random  vector  with  a  normal  distribution  of  which 
the  3<t  value  was  assumed  to  be  about  2  percent. 

Figs.  1-3  show  the  time  histories  of  the  attitude  com¬ 
mands  resolved  in  a  launch  inertial  frame  which  is  gener¬ 
ated  by  the  guidance  laws  for  the  simulated  flight  with 
the  nominal  condition.  Since  the  commands  generated  by 
‘PEG’,  ‘PEG  +  M.EG.’,  and  ‘PEG  +  L.D.C.’  are  very  sim¬ 
ilar,  those  generated  by  ‘PEG  +  M.I.G.’  are  omitted  in 
Figs.  1-3.  It  is  noticed  that  the  attitude  commands  gener¬ 
ated  by  different  guidance  schemes  show  little  difference 
for  the  second  stage  flight  with  a  short  duration.  However, 
those  generated  during  the  longer,  third  stage  flight  show 
some  differences  between  guidance  schemes.  For  the  yaw 
command,  ‘IGM’  appears  to  be  more  close  to  ‘Optimal’ 
than  ‘IGM  +  H.O.T’,  and  vice  versa  for  the  pitch  com¬ 
mand.  Also  as  can  be  expected,  ‘PEG’  and  ‘PEG  +  L.D.C.’ 
show  no  virtual  difference  and  other  states  omitted  here 
show  similar  results. 

Fig.  4  shows  that  the  time-to-go  prediction  error  of 
‘IGM’  gets  noticeably  larger  than  those  of  other  guidance 
schemes  with  a  longer  time-to-go.  It  also  shows  that  if  we 
consider  the  higher-order  terms,  this  error  can  be  reduced 
to  a  similar  level  with  that  of  PEG.  As  previously  men¬ 
tioned,  we  performed  five-hundred  different  simulations 
of  off-nominal  conditions  to  compare  the  guidance  schemes 
and  to  test  the  effect  of  inclusion  of  the  higher-order  terms 
to  the  original  IGM.  And  as  depicted  in  Fig.  4,  regardless 


of  the  scheme,  the  off-nominal  histories  represented  by  the 
500  thin  lines  show  similar  characteristics  as  the  nominal 
profile  represented  by  the  bold  line. 

Orbit  injection  accuracy  for  the  nominal  condition  with 
different  guidance  stop  conditions  is  shown  in  Table  2, 
where  the  orbit  injection  errors  are  defined  as  the  difference 
between  the  final  injection  states  obtained  from  the  simula¬ 
tions  and  the  target  orbit  values  used  for  guidance.  The 
guidance  laws  are  set  to  be  terminated  if  the  velocity-to- 
go  is  less  than  the  specified  value,  and  it  is  noticed  that 
‘IGM’  can  produce  some  erroneous  engine  shutdown  time 
with  an  early  termination  of  the  guidance  algorithm,  which 
again  may  result  in  orbit  injection  errors.  Attitude  com¬ 
mands  after  guidance  termination  are  usually  generated 
to  be  constant  or  to  change  linearly  with  time  and  here 
we  chose  the  linear  commands,  so  that  the  attitude  rates 
become  constant.  Fig.  5  illustrates  the  targeting  results  with 
the  off-nominal  conditions.  As  expected,  ‘IGM’  again  pro¬ 
duces  lower  performance  than  others  with  an  earlier  guid¬ 
ance  stop  and  ‘IGM  +  H.O.T.’  shows  similar  performance 
with  those  of  ‘PEG’s.  Fig.  6  shows  that  the  linear  angle 
assumption  of  IGM  in  Eqs.  (1)  and  (2)  is  more  valid  in 
‘IGM  +  H.O.T’  during  the  third  stage  flight  with  both 
the  off-nominal  and  nominal  conditions  although  there 
are  no  differences  between  them  in  the  short  flight  duration 
of  the  stage  2.  In  Fig.  7  it  is  observed  that  the  effects  of 
higher-order  nonlinear  terms  increase  with  remaining  flight 
time  for  the  simulated  trajectories. 

Table  2  and  Fig.  5  also  show  that  ‘PEG’,  ‘PEG  + 
M.I.G.’,  and  ‘PEG  +  L.D.C.’  have  little  differences  in  the 
orbit  injection  accuracy.  In  spite  of  the  similarity  in  the 
converged  guidance  commands  illustrated  in  Figs.  1-3, 
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Fig.  7.  Higher-order  terms  of  (a)  velocity  and  (b)  position  in  IGM  with  the  nominal  and  off-nominal  conditions. 


Fig.  8  illustrating  the  intermediate  values  during  the  itera¬ 
tion  process  shows  that  the  effect  of  inaccurate  initial  guess 
for  the  predictor-corrector  iteration  process  of  ‘PEG’ 
increases  as  the  time-to-go  approaches  0.  At  every  start 
of  the  iteration,  the  effect  of  incorrect  guesses  appears  as 
the  peak  points  in  blue  and  black  lines  in  Fig.  8.  If  the 
time-to-go  is  not  close  to  0,  after  a  small  number  of  itera¬ 
tions  they  converge  to  similar  values  with  those  of 
‘PEG  +  M.I.G.’  and  ‘PEG  +  L.D.C.’,  which  do  not  have 
the  initial  guess  effect.  However,  inaccuracy  in  the  initial 


guess  of  vgo  increases  as  the  time-to-go  reduces  to  0  and 
produces  incorrect  guidance  parameters  different  from 
those  of  ‘PEG  +  M.I.G.’  and  ‘PEG  +  L.D.C.’.  Fig.  9 
explains  the  reason  why  the  absolute  error  in  the  initial 
guess  of  vgo  assumed  by  the  magnitude  of  Avgra„  in  Eq. 
(45)  changes  slightly  during  the  flight  time,  while  its  relative 
error  increases  since  vgo  reduces  to  0  as  the  launch  vehicle 
approaches  the  target  point.  The  effect  of  incorrect  starting 
point  also  increases  with  the  duration  of  the  guidance  com¬ 
putation  cycle,  as  expected  from  Fig.  9.  An  earlier 
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Fig.  8.  Change  in  the  guidance  parameter  of  lx  in  PEGs  during  the  iteration  process  with  the  nominal  and  off-nominal  conditions.  This  parameter  is 
sensitive  to  the  initial  guess  of  ‘PEG’  and  the  sensitivity  also  increases  with  the  guidance  major  cycle  update  time  and  as  tgo  approaches  0. 


Fig.  9.  Initial  guess  error  relative  to  vgo.  The  relative  error  increases  as  tgo  approaches  0. 


divergence  is  observed  in  the  update  cycle  of  2.0  (sec)  when 
compared  with  the  0.5  (sec)  case  in  Fig.  8,  and  thus  an  early 
termination  of  the  guidance  is  required  to  prevent  the  insta¬ 
bility.  It  also  affects  the  number  of  the  iterations  to  be  per¬ 
formed  at  each  guidance  cycle,  as  illustrated  in  Fig.  10, 
from  which  it  is  observed  that  ‘PEG’  requires  more  itera¬ 
tions  than  other  algorithms  do  in  most  cases.  Here  the  iter¬ 
ation  is  set  to  be  stopped  when  the  difference  between  the 
updated  and  original  tgo  is  less  than  0.0001  (sec),  and  this 


criterion  can  be  mitigated  to  speed  up  the  computation. 
Figs.  8  and  10  also  include  off-nominal  results  obtained 
for  +3ff  perturbed  thrust  conditions  of  the  stage  2,  and 
they  also  show  similar  tendencies  with  the  nominal  condi¬ 
tions.  Also  we  checked  results  for  the  randomly  selected 
500  perturbed  conditions  and  observed  the  same  tenden¬ 
cies,  although  detailed  results  are  omitted  here. 

The  simulation  results  of  PEGs  imply  that  although  the 
original  predictor-corrector  algorithm  has  an  advantage 
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Time-to-go  (sec) 


Fig.  10.  Number  of  iterations  in  PEGs  with  the  nominal  and  off-nominal  conditions.  ‘PEG’  with  a  larger  initial  guess  error  results  in  more  iterations  than 
those  of  others. 


having  less  computational  complexity  than  the  differential 
corrector  algorithm  does  in  that  it  does  not  require  the 
derivatives  of  the  guidance  parameters,  the  convergence 
can  be  often  trickier  than  that  of  the  differential  corrector 
algorithm.  It  is  also  shown  that  a  slight  modification  of 
the  initial  guess  can  improve  the  convergence  characteris¬ 
tics  of  PEG. 

5.  Conclusions 

This  paper  considers  two  renowned  explicit  guidance 
algorithms  for  space  launch  vehicles,  IGM  and  PEG 
developed  in  the  space  era  of  1960-1 970s  and  presents 
performance  comparison  and  in-depth  analysis  of  each 
algorithm.  We  showed  that,  by  including  originally 
ignored  higher-order  nonlinear  terms  in  IGM,  the  perfor¬ 
mance  of  IGM  can  be  comparable  to  that  of  PEG.  And 
the  assumption  employed  in  the  derivation  of  IGM  that 
the  suboptimal  commands  are  in  linear  forms  becomes 
more  valid  with  modified  IGM.  Simulation  results  for 
both  nominal  and  off-nominal  conditions  imply  that 
accurate  representation  of  thrust  integrals  with  higher- 
order  nonlinear  terms  is  necessary,  especially  for  missions 
with  long  flight  time  and/or  range.  For  PEG,  we  showed 
that  sensitivity  to  the  incorrectly  guessed  velocity-to-go 
can  be  reduced  by  improving  initial  guess  for  the  itera¬ 
tive  predictor-corrector  computation  process.  Compari¬ 
son  with  the  linear  differential  corrector  approach  is 
performed  to  verify  the  proposed  method.  In  future 
study,  dependence  of  the  convergence  characteristics  of 
the  iteration  process  on  time-to-go  will  be  analyzed  in 
detail. 
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